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Pj ■ Abstract 
cn . 

We summarize the status of next-to-leading order perturbative quantum 
chromodynamics (pQCD) calculations of the cross sections for the produc- 
tion of squarks, gluinos, neutralinos, charginos, and sleptons as a function of 
the produced sparticle masses in proton-antiproton collisions at the hadronic 



center-of-mass energy 2 TeV. 



I. PERTURBATIVE QCD RESULTS 

The possibility of supersymmetry (SUSY) at the electroweak scale and the ongoing search 
for the Standard Model (SM) Higgs boson constitute two major related aspects of the 
motivation for the Tevatron upgrade currently under construction at Fermilab. The increase 
in the center-of-mass energy to 2 TeV and the luminosity to an expected 2 fb _1 , together with 
detector improvements, should permit discovery or exclusion of supersymmetric partners of 
the standard model particles up to much higher masses than at present 

Estimates of the production cross sections for pairs of supersymmetric particles may 
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be computed analytically from fixed-order quantum chromodynamics (QCD) perturbation 
theory Calculations that include contributions through next-to-leading order (NLO) in 
QCD have been performed for the production of squarks and gluinos 0, top squark pairs 
||, slepton pairs ||, gaugino pairs 0, and the associated production of gauginos and gluinos 
H. The cross sections can be calculated as functions of the sparticle masses and mixing 
parameters. 

In a recent paper [f7|, Berger, Klasen, and Tait provide numerical predictions at next-to- 
leading order for the production of squark-antisquark, squark-squark, gluino-gluino, squark- 
gluino, and top squark - antitop squark pairs in proton-antiproton collisions at the hadronic 
center-of-mass energy 2 TeV. These calculations are based on the analysis of Refs. , and 
the CTEQ4M parametrization || of parton densities. The hard scale dependence of the 
cross section at leading order (LO) in perturbative QCD is reduced at NLO but not absent. 
An estimate of the theoretical uncertainty at NLO is approximately ±15 % about a central 
value. The central value is obtained with the hard scale chosen to be equal to the average 
of the masses of the produced sparticles, and the band of uncertainty is determined from a 
variation of the hard scale from half to twice this average mass. The next-to-leading order 
contributions increase the production cross sections by 50 % and more from their LO values. 
For example, in the case of squark-antisquark production the next-to-leading order cross 
section lies above the leading order cross section by 59 %. This increase translates into a 
shift in the lower limit of the produced squark mass of 19 GeV. The cross sections for squark- 
antisquark production, gluino pair production, and the associated production of squarks and 
gluinos of equal mass are of similar magnitude, whereas the squark pair production and top 
squark-antitop squark production cross sections are smaller by about an order of magnitude 

HI- 

The cross sections reported in Ref. |7[ are for inclusive yields, integrated over all trans- 
verse momenta and rapidities. In the search for supersymmetric states, a selection on trans- 
verse momentum will normally be applied in order to improve the signal to background 
conditions. The theoretical analysis can also be done with similar selections. A tabulation 
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of cross sections for various squark and gluino masses is available upon request from the 
authors of Ref . |7|] . 

Next-to- leading order calculations of the production of neutralino pairs, chargino pairs, 
and neutralino-chargino pairs are reported to be on the way to completion ||, but final 
numerical predictions are not yet available for general use. 

The strongly interacting squarks and gluinos may also be produced singly in association 
with charginos and neutralinos. Leading-order production cross sections for the associated 
production of a chargino plus a squark or gluino and of a neutralino plus a squark or 
gluino are published || , and a next-to- leading order calculation of associated production of 
a gaugino plus a gluino is now available |J. 

Berger, Klasen, and Tait || compute total cross sections for all the gaugino-gluino pro- 
duction reactions gx\i-4) an d gXn-2) m next-to-leading order SUSY-QCD. For numerical 
results, they select an illustrative mSUGRA scheme in which the GUT scale common scalar 
mass m = 100 GeV, the common gaugino mass m 1( / 2 = 150 GeV, the trilinear coupling 
A = 300 GeV, tan(/3) = 4 and sgn(yu) = +. (The sign convention for A is opposite to 
that in the ISASUGRA code). They convolute the NLO hard partonic cross sections with 
the CTEQ4M parametrization || of parton densities, and present physical cross sections 
as a function of the g mass or of the average mass m = (m^ + rrig)/2. For pp collisions at 
\/~S = 2 TeV the cross sections at rrig = 300 GeV range from (9(lpb) for the x\ anc ^ the 
Xf to C(10~ 3 pb) for the The gx°{i 2) an d gxt cross sections are of hadronic size despite 
the fact that the overall coupling strength is 0(otEw a s) not 0(a 2 s ). The masses of the Xn 2) 
and gxf are significantly smaller in a typical mSUGRA scenario than those of the squarks 
and gluinos. The phase space and the parton luminosity are therefore greater for associated 
production of a gluino and a gaugino than for a pair of squarks or gluinos, and the smaller 
coupling strength is compensated. The next-to-leading-order cross sections are enhanced 
by typically 10% to 25% relative to the leading order values. The theoretical uncertainty 
resulting from variations of the factorization/renormalization scale is approximately ±10% 
at NLO for the x% an d the Xii a fector of 2 smaller than the LO variation. Shown in Fig. 1 
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are the predicted cross sections as a function of the average mass. 

Baer, Harris, and Reno Q compute total cross sections for all the slepton pair production 
reactions Sl^l, &r&r and Vl&l in next-to-leading order QCD. The analytic calcula- 

tions are very similar to the QCD corrections to the Standard Model massive lepton-pair 
production (Drell-Yan) process. Numerical results are based on the CTEQ4M parametriza- 
tion H of parton densities. For pp collisions at = 2 TeV, the cross sections range from 
O(lpb) at m s i ep ton = 50 GeV to O(10 _3 pb) at m s i cpto n = 200 GeV. The next-to-leading-order 
cross sections are enhanced by typically 35% to 40% relative to the leading order values. 
The theoretical uncertainty resulting from variations in the hard scattering scale and par- 
ton distribution functions is approximately ±15%. In the mSUGRA model, slepton pair 
production is most important for small values of the parameter m . The next-to-leading 
order enhancements of slepton pair cross sections at Tevatron energies can push predictions 
for leptonic SUSY signals to higher values than typically quoted in the literature in these 
regions of model parameter space. 

For current expectations of the hierarchy of masses and cross sections, consult Ref. ||T|. 



II. MONTE CARLO METHODS 

Experimental searches for supersymmetry rely heavily on Monte Carlo simulations of 
cross sections and event topologies. Two Monte Carlo generators in common use for hadron- 



hadron collisions include SUSY processes; they are IS A JET |10[ and SPYTHIA |TT],|T2[. Both 
the Monte Carlo approach and the fixed order pQCD approach have different advantages and 
limitations. Next-to-leading order perturbative calculations depend on very few parameters, 
e.g., the renormalization and factorization scales, and the dependence of the production cross 
sections on these parameters is reduced significantly in NLO with respect to LO. Therefore, 
the normalization of the cross section can be calculated quite reliably if one includes the NLO 
contributions. On the other hand, the existing next-to-leading order calculations provide 
predictions only for fully inclusive quantities, e.g., a differential cross section for production 
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of a squark or a gluino, after integration over all other particles and variables in the final 
state. In addition, they do not include sparticle decays. This approach does not allow for 
event shape studies nor for experimental selections on missing energy or other variables 
associated with the produced sparticles or their decay products that are crucial if one wants 
to enhance the SUSY signal in the face of substantial backgrounds from Standard Model 
processes. 

The natural strength of Monte Carlo simulations consists in the fact that they generate 
event configurations that resemble those observed in experimental detectors. Through their 
parton showers, these generators include, in the collinear approximation, contributions from 
all orders of perturbation theory. In addition, they incorporate phenomenological hadroniza- 
tion models, a simulation of particle decays, the possibility to implement experimental cuts, 
and event analysis tools. However, the hard-scattering matrix elements in these generators 
are accurate only to leading order in QCD, and, owing to the rather complex nature of 
infra-red singularity cancellation in higher orders of perturbation theory, it remains a diffi- 
cult challenge to incorporate the full structure of NLO contributions successfully in Monte 
Carlo simulations. The limitation to leading-order hard-scattering matrix elements leads 
to large uncertainties in the normalization of the cross section. The parton shower and 
hadronization models rely on tunable parameters, another source of uncertainties. 

In Ref. a method is suggested to improve the accuracy of the normalization of cross 
sections computed through Monte Carlo simulations. In this approach, the renormalization 
and factorization (hard) scale in the Monte Carlo LO calculation is chosen in such a way 
that the normalization of the Monte Carlo LO calculation agrees with that of the NLO 
perturbative calculation. The scale choice depends on which partonic subprocess one is 
considering and on the kinematics. This choice of the hard scale will affect both the hard 
matrix element and the initial-state and final-state parton shower radiation. On the other 
hand, an alternative rescaling of the cross section by an overall i^-factor will have no bearing 
on the parton shower radiation. A reduction in the hard scale leads generally to less evolution 
and less QCD radiation, and vice-versa, in the initial- and final-state showering. A change 
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of the hard scale will be reflected in the normalization of the cross section as well as in the 
event shape. Investigations are underway to determine how significant the changes are in 
computed final state momentum distributions. 



6 



REFERENCES 



[1] M. Carena, R.L. Culbertson, S. Eno, H.J. Frisch, and S. Mrenna, "The Search for 
Supersymmetry at the Tevatron Collider", Argonne report ANL-HEP-PR-97-98, [hep- 
ex/9712022t to be published in Rev. Mod. Phys. 



[2] W. Beenakker, R. Hopker, M. Spira, and P.M. Zerwas, Nucl. Phys. B492, 51 (1997). 

[3] W. Beenakker, M. Kramer, T. Plehn, M. Spira, and P.M. Zerwas, Nucl. Phys. B515, 3 
(1998). 

[4] H. Baer, B.W. Harris, and M.H. Reno, Phys. Rev. D 57 (1998) 5871. 

[5] W. Beenakker, M. Klasen, M. Kramer, T. Plehn, M. Spira, and P.M. Zerwas, Talks 
given by M. Kramer and M. Spira at the Workshop on Theory of LHC Processes, 
CERN, Geneva, Switzerland, February 1998. 

[6] E. L. Berger, M. Klasen, and T. Tait, Argonne report ANL-HEP-PR-99-03, 
|hep-ph/9gD235q . 

[7] E. L. Berger, M. Klasen, and T. Tait, Argonne report ANL-HEP-PR-98-48, 
|hep-ph/9807230| , to be published in Phys. Rev. D. 

[8] H.L. Lai, J. Huston, S. Kuhlmann, F. Olness, J. Owens, D. Soper, W.K. Tung, and H. 
Weerts, Phys. Rev. D 55, 1280 (1997). 

[9] S. Dawson, E. Eichten, and C. Quigg, Phys. Rev. D 31, 1581 (1985); H. Baer, 
D. D. Karatas, and X. Tata, Phys. Rev. D 42, 2259 (1990). 

[10] F.E. Paige, S.D. Protopopescu, H. Baer, and X. Tata, Brookhaven report BNL-HET- 
98-18, |hep-ph/980432Tl . 



[11] T. Sjostrand, Comput. Phys. Commun. 82 , 74 (1994). 
[12] S. Mrenna, Comput. Phys. Commun. 101,232 (1997). 



7 



FIGURES 
pp -> g% at VS = 2.0 TeV 
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FIG. 1. Total hadronic cross sections for the associated production of gluinos and gauginos at 
Run II of the Tevatron from Ref.[6]. NLO results are shown as solid curves, and LO results as 
dashed curves. The chargino cross sections are summed over positive and negative chargino states. 
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